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Adiabatic Perturbations in Homologous Conventional 
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ABSTRACT 

We perform a non-radial adiabatic perturbation analysis on homologous conventional 
polytropic stellar core collapses. The core collapse features a polytropic exponent T = 
4/3 relativistic gas under self-gravity of spherical symmetry while three-dimensional 
perturbations involve an adiabatic exponent 7 with 7 / T such that the Brunt- Vaisala 
buoyancy frequency Af does not vanish. With proper boundary conditions, we derive 
eigenvalues and eigenfunctions for different modes of oscillations. In reference to stellar 
oscillations and earlier results, we examine behaviours of different modes and the 
criterion for instabilities. The acoustic p— modes and surface f— modes remain stable. 
For 7 < r, convective instabilities appear as unstable internal gravity g~— modes. For 
7 > r, sufficiently low-order internal gravity g + — modes are stable, whereas sufficiently 
high-order g + — modes, which would have been stable in a static star, become unstable 
during self-similar core collapses. For supernova explosions, physical consequences of 
such inevitable g— mode instabilities are speculated. 
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1 INTRODUCTION 

Stability properties of core collapses in massive progenitor 
stars before supernova (SN) explosions (Goldreich & Weber 
1980 - GW hereafter; Goldreich et al. 1996; Lai 2000; Lai 
& Goldreich 2000; Blondin et al. 2003; Murphy et al. 2004; 
Burrows et al. 2006, 2007; Lou & Cao 2008; Cao & Lou 
2009) have come to focus after three decades, because of the 
realization after numerous unsuccessful one-dimensional SN 
simulations that the breakdown of spherical symmetry in- 
evitably occurs and plays a key role in SN explosions (e.g. 
Burrows 2000). Numerical simulations (e.g. Bruenn 1989a, 
b) indicate that pre-SN stellar core collapses may be approx- 
imately described by a homologous process, first analyzed by 
GW for a conventional polytropic equation of state (EoS) 
P — Kp F where P and p are the pressure and density and 
both k, and F = 4/3 are constant. For perturbations obeying 
identical EoS of the background core collapse, GW explored 
linear stability properties of such stellar collapses and con- 
cluded that these collapses are stable. Yahil (1983) extended 
homologous collapses to polytropic exponent F < 4/3, not- 
ing the presence of an outer supersonic envelope besides the 
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inner core collapse. Lai (2000) performed perturbation anal- 
ysis to these extended solutions and claimed that perturba- 
tions are stable for T > 1.09. 

Lou & Cao (2008) substantially extended the T = 4/3 
self-similar solutions, including those of GW, using a gen- 
eral polytropic EoS with a temporally and radially variable 
n being conserved along streamlines. We obtained a broad 
family of homologous core collapses allowing n as an arbi- 
trary function of the independent self-similar variable. 

The specific entropy of an ideal or perfect gas is 



(7-1) 



In 



+ constant , 



(1) 



where ks is the Boltzmann constant and 7 is the adia- 
batic exponent with 7 / T. Here F determines the struc- 
ture and evolution of core collapses while 7 controls specific 
entropy perturbation properties. The earlier isentropic as- 
sumption requires 7 = V (e.g. GW). Numerical simulations 
(e.g. Bethe et al. 1979; Bruenn 1985, 1989a, 1989b; Woosley 
et al. 1993, 2002) for the structure of massive stars and SN 
explosions support variable radial distributions of specific 
entropy. Thus the model of Lou & Cao (2008) allows an 
arbitrary radial profile of specific entropy. Using this gen- 
eral polytropic model, we conducted a non-radial adiabatic 
perturbation analysis and classified different perturbation 
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modes parallel to stellar oscillations (see Unno et al. 1979 
and Cao & Lou 2009) and concluded that in addition to 
internal gravity g~— modes for convective instabilities, suf- 
ficiently high-order g + — modes also become unstable. 

Cao & Lou (2009) demonstrated that the criterion of 
using the sign of the Brunt- Vaisala, buoyancy frequency 
squared AT 2 (see definition [I5j for the existence of g— modes 
remains valid in self-similar dynamic core collapses. In a con- 
ventional polytropic process, isentropic perturbations with 
7 = r make all g— modes disappear, i.e. Af 2 = 0. Thus previ- 
ous perturbation analyses (GW; Lai 2000; Lai & Goldreich 
2000) considered only the stability of acoustic modes. In ad- 
dition to p— modes and f— modes, Cao & Lou (2009) found 
g— mode instabilities in stellar core collapses and speculated 
possible consequences for neutron star kicks etc. 

From the existence criterion of Af 2 ^ for g— modes, 
we realize that in a conventional polytropic core collapse, 
non-radial adiabatic perturbations with 7 / T should allow 
g— modes. Such nonisentropic assumption has been invoked 
decades ago (e.g. Ledoux 1965; Unno et al. 1979) and ap- 
plied to helioseismology and white dwarf oscillations (e.g. 
Shibahashi et al. 1988). For star formation, McKee & Holli- 
man (1999) also studied nonisentropic radial perturbations 
in molecular clouds. They used a thermal free energy to dis- 
cuss cloud stability and gave critical cloud masses for locally 
and globally adiabatic models. These early results motivate 
us to perform such a perturbation analysis for stellar core 
collapses and study the existence and stability of g— modes 
in addition to other modes. The main thrust of this Letter 
is to show such g— mode instabilities and speculate conse- 
quences for stellar core collapses and SN explosions. 



where the free-fall timescale tff is defined by 
1 -1/2 



tff = 



(7) 



and p is a dimensionless 'frequency'. We do not use stream 
function for velocity perturbations to allow for vorticities. 
The dynamic collapse background is that of GW. Defining 
a dimensionless collapse parameter A (Yahil 1983) as 



1 = SirGpar 2 
A ~ 
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(8) 



where the subscript indicates background variables and 
Uff is the free-fall velocity, we obtain the background density 
profile determined by the equation for f(x), namely 



_i_d_ 

x 2 dx 



>df(x) 



dx 



+ f = A , 



(9) 



with boundary conditions /(0) = 1 and /'(0) = 0. Other 
background variables are then readily determined. Solutions 
of eq (0 are sensible for core collapses of massive stars 
when density vanishes at a moving boundary Xb- For this 
purpose, the variable a should be negative and the time re- 
versal operation must be taken in the governing equations. 
We assume adiabatic perturbations with the exponent 7 dif- 
ferent from r = 4/3. In our model, we define vi = pwi 
and m = p[p + (A/2) 1//2 ] for convenient analysis. After the 
linearization of hydrodynamic equations including conser- 
vations of mass and momentum, Poisson equation for self- 
gravity and the adiabatic EoS, the angular variation factors 
of wi, (3i, fi and ipi can be separated out involving the 
spherical harmonics Yj m (#, (f>). We further let wi take the 
consistent form of 



2 PERTURBATIONS IN CORE COLLAPSES 

The governing equations include conservations of mass and 
momentum, Poisson equation for self-gravity and EoS. The 
gas is ideal and the background is a V — 4/3 conventional 
polytrope while the adiabatic perturbation with 7 7^ T 
conserves the specific entropy along streamlines. The back- 
ground collapse is spherically symmetric and as in GW, we 
introduce a dimensionless independent variable x = r/a(t) 
where a(t) = p c 1 ^ 3 [k/(tvG)] 1 ^ 2 is the Jeans length with a 
central mass density p c (t) and G is the gravitational con- 
stant. Physical variables are cast into the following forms 
where the first term of each variable is the dynamic back- 
ground and the last term is the first-order perturbation: 



• / \ a (t) , 
u = a(t)x H —vi(x, 

tff 
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Here the time-dependent temporal factor r(t) is 



r{t) = exp (p J t f jdt'\ , 



(4) 
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wi = e r w r (x,t)Y lm +w t (x,t)Vj_Yi m + V± x (w rot Yi m e r ) , 

(10) 

where the modified transverse gradient operator is 



- d Id 

00 sm aq> 



(11) 



As Wrot decouples from other perturbation equations, we 
suppress it in the following analysis. Eliminating f3i and fi, 
we obtain a set of ordinary differential equations (ODEs): 

/ d ( 2 , l(l + l)f 4 df 

[X Wr) Wt H Wr 



x 2 dx 



1 d ( 2 #l\ 1(1 + 1) 



7 dx 
(mxwt + ipi) = , (12) 
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dx 



— / (mxw t + ipi) , (13) 
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d . . (37 - 4) df , 
mw r — m -j^ \xw t ) = — — — — (mxwt + tpi) 



+ 



3(4 -3 7 ) (df 



if \dx 



w r , (14) 



where the same notation is used for the a;— dependent factor 
of ip\, ODEs (|12[l — (|14|) for non-radial perturbations reduce 
to those of GW by setting 7 = F = 4/3 and introducing a 
stream function for velocity perturbations. 
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Regular boundary conditions are imposed, namely 



if} i oc x , 



ipi OC X 



-O+i) 



w r — lw t 



3w r df /dx — mxwt 



for x — + + 



for x 



The last of the outer boundary conditions requires a zero 
Lagrangian pressure. With these boundary conditions, we 
can prove the orthogonality of eigenfunctions, i.e. J / 3 w^' • 
w'Vr = where superscripts and ' J ' indicate eigen- 
functions of different eigenvalues. Meanwhile, the variational 
principle (Chandrasekhar 1964) can be formulated for this 
eigenvalue problem. 

To solve this eigenvalue problem, we use the following 
numerical procedure. After solving the background physi- 
cal variables by using the fourth-order explicit Runge-Kutta 
scheme, we discretize ODEs (1 1 2 p — (|14l) with a proper mesh. 
Then the inverse iteration method (e.g. Wilkinson 1965) is 
applied to determine eigenvalues and eigenfunctions. To en- 
hance the efficiency of the inverse iteration method, a relax- 
ation scheme is also implemented. 



3 RESULTS OF OUR MODEL ANALYSIS 

Wc first test our numerical code to solve 7 = 4/3 isentropic 
cases and confirm the results of Cao & Lou (2009) to guar- 
antee the correctness and reliability of the code. This check 
also confirms the errors of p— mode eigenvalues in GW. 

The spectra of eigenvalues versus the mode degree I 
contain characteristic features of oscillations. A typical spec- 
trum of 7 > 4/3 modes is shown in Fig. [1] In these cases, 
all eigenvalues m are negative. We use two dashed lines 
to separate three distinct classes of modes: p— modes for 
I > 0, a unique branch of f— modes for I > 2 and g + — modes 
for / > 1 (see Cowling 1941 for mode classifications of 
stellar oscillations and Unno et al. 1979 for properties of 
different modes.). The p— modes are pressure-driven and 
gravity-modified acoustic oscillations. The f-modes are es- 
sentially trapped surface Lamb waves (Lamb 1932; Unno et 
al. 1979; Lou 1990, 1991) and decay exponentially inwards. 
The g— modes are driven by buoyancy and are trapped deep 
inside the core. The g + — modes are one subclass of g— modes 
with 7 > 4/3. For a given I, the absolute value of eigenvalue 
\m\ decreases as the radial order increases. In that limit, 
eigenvalue m approaches 0. There is no I = g + — mode. 
The other subclass g~— modes exist for 7 < 4/3 with sim- 
ilar characteristics of g + — modes, except that the eigenval- 
ues m are positive and decrease towards zero as the radial 
order increases. By numerical explorations, we find similar 
behaviours of each perturbation mode in the core collaps- 
ing background as those in stellar oscillations. As p— modes 
have been studied by GW for 7 = F = 4/3, we emphatically 
show g + — mode eigenfunctions and their eigenvalues in Fig. 

m 

Between p— modes and g + — modes exists the unique 
branch of f— modes for / > 2. Eigenvalues of f— modes sepa- 
rate the p— and g— modes. Such f— modes are characterized 
by eigenfunctions of density and radial velocity perturba- 
tions without nodes. They are acoustic in nature and relate 
to the surface Lamb waves (Lamb 1932; Lou 1990, 1991). 

In stellar oscillations, the existence criterion for both 
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Figure 1. Sample spectra of eigenvalues m versus spherical har- 
monic degree I are shown with the parameter A = 0.002 and the 
adiabatic exponent of perturbations 7 = 1.6. All eigenvalues lie in 
the m < regime. Two dotted lines divide the region into three 
distinct zones for eigenvalues of different modes. There is a sign 
difference between the definitions of the eigenvalue m here and in 
stellar oscillations (see figures of Cox 1976 for a comparison). 
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Figure 2. Eigenfunctions of g + — modes with 1,2,3 nodes from 
bottom to top panels with eigenvalues m = —0.0249, —0.0135, 
—0.0085, respectively, are presented here. The solid and dashed 
curves represent xf 5 ^w r and x 2 f s / 4 wt, respectively. Factors 
a;/ 5 / 4 and x 2 / 5 / 4 are multiplied for compact presentations. 



types of g— modes depends on the square of the Brunt- 
Vaisala buoyancy frequency Af 2 defined by 



A/" 2 = — ~4— 
p dr 



dlnp 

dr 



ldlnP 

7 



dr 



(15) 



where P(r) and p(r) are for a hydrostatic equilibrium, and 7 
is the adiabatic exponent of perturbation. Since our numeri- 
cal explorations reveal that if an eigenvalue m with its eigen- 
function can be obtained for A = for the hydrostatic limit, 
its counterpart can also be determined for A > in a con- 
tinuous manner. Thus, the existence criterion for g— modes 
in stellar oscillations can be applied in our analysis by using 
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Figure 3. Displayed are eigenvalue m variations of the three 
lowest acoustic p— modes and internal gravity g + — modes of I = 
1, 2 and of surface f— modes for I = 2 with the collapse parameter 
A characterizing different self-similar dynamic backgrounds. The 
adiabatic exponent of perturbations is 7 = 1.5. 



Figure 4. Displayed are eigenvalue m variations of the two lowest 
acoustic p— modes and internal gravity g + — modes and of surface 
f— modes for I = 2 with the adiabatic exponent 7 varying from 
1.34 to 1.66. The collapse parameter is A = 0.002 for the conven- 
tional polytropic stellar core collapse of T = 4/3. 



partial derivatives instead of derivatives with respect to the 
radius. For 7 < 4/3 (TV 2 < 0), g~ — modes appear, whereas 
for 7 > 4/3 (TV 2 > 0), g + -modes manifest. For 7 = T = 4/3 
(A/" 2 = 0), both g— modes are suppressed. 

We have readily computed the p— mode eigenfunctions 
with 7 / T = 4/3 and they appear qualitatively similar to 
those shown in GW and in Cao & Lou (2009). 

We demonstrate mathematically that a g + — mode does 
not change to a g~ — mode or vice versa as A increases from 
to Am, i.e. for 7 7^ 4/3, there is no m — eigenvalue for 
all < A < Am- By setting m = in ODEs {T2j— (dU) and 
after straightforward manipulations, we derive 

/ [|V(^iy im )| 2 +3/ 2 ^Fi 2 m ]rf 3 r = , (16) 
Jv 

which requires tpi = 0. Thus no non-trivial eigenfunctions 
can be found. Consequently, the eigenvalue curve m(A) for 
each perturbation mode does not intersect the m = line. 
We show the variation trends of eigenvalues m as A increases 
from to Am in Fig. [3] for typical eigenvalues m of the three 
lowest p— modes and g— modes for I = 1, 2 and of f— mode 
for / = 2 as functions of A. Fig. [3] shows that eigenvalues m 
larger than 25 Aa/ /8 approach this limiting value as A — » \m 
(see GW). Fig. [3] illustrates the variation of eigenvalues m 
for the first two lowest orders of p— modes and g + — modes 
as well as that of f— modes all for A = 0.002 and I = 2 as the 
adiabatic exponent 7 varies from 1.34 to 1.66. 

A key question for perturbation analysis is the stability 
of a core collapse. The inequality TV 2 < for stellar pertur- 
bations is the Schwarzschild criterion for convection. In our 
dynamic core collapse background, we emphasize that this 
inequality is only sufficient but not necessary for instability. 
The temporal factor r(t) for A > bears a power law F 
with £ = -1/6 ± [1/36 + 2m/(9A)] 1/2 . For a time reversal 
operation, GW showed that the compression of a collapse is 
responsible for an amplification of t -1 ^ 6 . Thus the determi- 
nant under the square root in £ decides the mode stability. 
For m > —A/8 so that £ has two real roots, the mode is 



unstable for one of the two roots. For m < —A/8 so that £ 
has a pair of complex conjugate roots, a perturbation oscil- 
lates stably in addition to the compression amplification of 
core collapse. By this new criterion for instability, we find 
that all g~— modes and sufficiently high-order g + — modes 
are unstable. It can be shown by a local analysis that the 
definition of J\f remains valid for dynamical core collapses, 
such that g~— modes lead to convective instabilities. It also 
announces that g + — mode instabilities are uniquely associ- 
ated with self-similar dynamic core collapse. 

We compare our results with earlier analyses (GW; Lai 
2000; Cao & Lou 2009). First, in contrast to oscillations in 
a static star, we examine non-radial adiabatic perturbations 
in a self-similar conventional polytropic core collapse. The 
adiabatic exponent 7 differs from the background F = 4/3, 
whereas earlier stability analyses of the dynamic background 
are restricted to 7 = F — 4/3. Secondly, for a conven- 
tional polytropic core collapse, we classify different adia- 
batic perturbation modes, including p— modes, g— modes 
and f— modes. GW and Lai (2000) studied only p— mode 
perturbations as their isentropic EoS makes g— modes van- 
ish. Cao & Lou (2009) obtained g— modes by a general poly- 
tropic EoS with a variable specific entropy distribution, even 
though the perturbation EoS remains the same as the back- 
ground EoS. In comparison, our model here describes non- 
radial adiabatic perturbations with 7 / T = 4/3. This leads 
to N 7^ so that g— modes may exist. We indeed confirm this 
by numerical explorations. For various efficiencies of heat 
transport and radiative losses, the case of 7 = V should be a 
special and rare situation. Therefore g— modes should exist 
in stellar core collapses in general. 

For core-collapse SN explosions, such g— mode instabil- 
ities should bear physical consequences. These instabilities 
complement those revealed by Cao & Lou (2009). They oc- 
cur during the pre-SN core collapse stage and should also 
influence the formation of proto-neutron stars and subse- 
quent emergence and evolution of rebound shocks. At least, 
such instabilities lead to early breakdown of spherical sym- 
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metry before the emergence of rebound shocks (e.g. Lou & 
Wang 2006, 2007; Wang & Lou 2008; Hu & Lou 2009). So 
far, most proposed instabilities occur either in a massive 
progenitor star prior to the core collapse (e.g. Goldreich et 
al. 1996) or after the core rebound (e.g. Blondin et al. 2003; 
Burrows et al. 2006, 2007). Such instabilities may affect the 
formation, motion and evolution of a proto-neutron star or 
a pulsar as speculated by Cao & Lou (2009). Moreover, our 
perturbation analysis serves to link perturbations in a mas- 
sive progenitor star and those after the core rebound. For 
example, the 'e-mechanism' (e.g. Goldreich et al. 1996) may 
provide seed perturbations during a core collapse while such 
perturbations during a core collapse may stimulate those 
after the core rebound. 

We emphasize that while p— modes, f— modes and low- 
order g + — modes are stable, if radiative losses and diffusive 
processes are involved instead of the adiabatic approxima- 
tion, these modes might become overstable. 



4 SUMMARY AND CONCLUSIONS 

In this Letter, we examine 3D adiabatic perturbations in 
a self-similar conventional polytropic collapsing core with 
7 7^ r = 4/3. The gas is non-isentropic in the sense that the 
adiabatic index 7 of perturbations differs from that of the 
background F = 4/3. As 7 determines the specific entropy 
perturbation conserved along streamlines, the non-isentropic 
process actually involves a nonzero buoyancy frequency jV 
giving rise to g— mode perturbations. In comparison, pertur- 
bation analysis of Cao & Lou (2009) emphasizes a general 
polytropic core collapse with a variable specific entropy (Lou 
& Cao 2008) and 3D adiabatic perturbations of 7 = F — 4/3. 

By imposing proper boundary conditions, we solve the 
perturbation eigenvalue problem and derive distinct modes: 
acoustic p— modes, surface f— modes, internal gravity g~ — 
and g + — modes which are classified by their series of eigen- 
values m and eigenf unctions in reference to stellar oscilla- 
tions. In parallel, g— modes involve two types, viz. g + — and 
g~ — modes; their existence depends on the square of the 
Brunt- Vaisala buoyancy frequency A/" 2 : the former requires 
TV 2 > while the latter needs JV' 2 < 0. For adiabatic pertur- 
bations, g + — modes correspond to 7 > 4/3 while g~— modes 
occur for 7 < 4/3. 

Stability properties of these modes in stellar core col- 
lapses are examined. The instability criterion shifts from 
m > in a static Lane-Emden polytropic sphere of A = 
to m > —A/8 for < A < Am core collapses. Consequently, 
p— modes, f— modes and sufficiently low-order g + — modes os- 
cillate stably. The g~ —modes are unstable leading to convec- 
tions. Sufficiently high-order g + — modes which would have 
been stable in a static polytrope now become unstable for 
< A < A a/ GW dynamic core collapse. The specific ra- 
dial order that g + — modes become unstable depends on the 
parameter pair of A and 7. 

We speculate that such inevitable g— mode instabilities 
may offer valuable clues to SN simulations and that the for- 
mation, motion and evolution of proto-neutron stars can be 
nontrivially influenced (Cao & Lou 2009). These perturba- 
tions again lead to instabilities before the core rebound and 
they serve as seeds of later fluctuations. Their excitations 
can come from oscillations of the progenitor star before the 



onset of core collapse. According to our numerical explo- 
ration and typical stellar parameters with 7 > Y = 4/3, 
unstable g + — modes appear with fairly high radial orders. 
This implies that the central mass blob can be quite small 
and a SN might even break the core into pieces without 
forming a NS. 
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